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It is argued that the d x 2_ y 2 superconductor exhibits the transition into d x 2_ y 2 + id xy state in 
the presence of magnetic impurities in certain conditions. Linear coupling between impurity mag- 
netization and the d xy component drives this transition. Violation of time-reversal symmetry and 
parity occurs spontaneously via the second order transition. In the ordered phase both impurity 
magnetization and d xy component of the order parameter develop and are proportional to each 
other. 
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It is well known that magnetic impurities destroy 
the singlet superconducting state due to spin scattering 
which breaks pair singlets Q . In the case of the gapless 
(with the nodes of the gap) d-wave superconductor both 
magnetic and nonmagnetic impurities produce a finite 
density of states at zero energy. These effects are a simple 
and direct consequence of the lifetime effects produced by 
impurities. These are well known "incoherent" effects of 
impurities in unconventional superconductors. After re- 
cent experiments by Movshovich et.al. we are led to 
believe that another phenomena is possible, namely the 
transition to the second superconducting phase as a re- 
sult of condensate interactions with magnetic impurities. 
The time reversal violating state is formed at low energy 
and the order parameter of the new phase is d x 2_ y 2 +id xy 
(d+id). In this phase the impurity spins acquire nonzero 
spin density along z-axis, i.e. out of plane. 

The physical origin of the instability comes from the 
fact that the d + id state has an orbital moment which 
couples to the magnetic impurity spins. In the pure phase 
one can think of d-wave state as an equal admixture of 
the orbital moment L z = ±2 pairs: 



Aq(0) = A o cos20 = ^(exp(2i<j 



+ exp(-2z<£)) (1) 



Here <j) is the 2D planar angle of the momentum on the 
Fermi surface, Ao is the magnitude of the d x 2_ y 2 compo- 
nent. We consider 2D d x 2_ y 2 superconductor, motivated 
by the layered structure of the cuprates. In the presence 
of the (ferromagnetically) ordered impurity spins S z the 
coefficients of the L z — ±2 components will shift linearly 
in S z with opposite signs: 

Ao(0)->^r((l+0S,)exp(2i0) 
+(l-sS,)exp(-2i0))=Ao(0)+t SMfl (2) 

where Ai(</>) = g/2 sin 20 - is the d xy component and g 
is the coupling constant. The relative phase 7r/2 of these 
two order parameters comes out naturally because d + id 
state has a noncompensated orbital moment L z = +2. 



This natural coupling of the orbital moment of the pairs 
to the spin magnetic moment leads to the instability of 
the d x 2_ y 2 state towards time reversal violating d + id 
state. 

Recent experimental observation of the surface- 
induced time-reversal violating state in YBCO suggests 
that the secondary component of the order parameter 
(d+is) can be induced jjj. Theoretical explanation, based 
on surface-induced Andreev states has been suggested by 
Sauls and co-workers ||]. The source of the secondary 
component is the bending of the original d x 2_ y 2 order 
parameter at the surface. 

In a different approach Laughlin || argued that the 
d x 2_ y 2 state is unstable towards d + id state in the bulk 
in the perpendicular magnetic field at low enough tem- 
peratures. The time reversal and parity are broken by 
external field in this case. The linear coupling of the sec- 
ondary order parameter to the external field is central 
to his consideration and results in the first order phase 
transition into d+id state. This transition was suggested 
to be responsible for the kink-like feature in the thermal 
conductivity in experiments by Krishana et.al. Q 

The purpose of this paper is to show that that there 
is yet another possibility to break time reversal and par- 
ity: spontaneously in the bulk of the d-wave state due 
to coupling to the impurity spins. We find that: 1) 
Original d x 2_ y 2 is unstable towards formation of the bulk 
d x 2_ y 2 +id xy phase in the presence of the magnetic impu- 
rities under certain conditions. Time reversal symmetry 
and parity (T and P) are spontaneously broken as the 
result of interaction between conduction electrons and 
impurity spins. Below we consider magnetic impurities 
with nonzero spin 5, such as Ni with S=l at low tem- 
peratures. In the ordered phase both the spontaneous 
magnetization of impurity spins and second component 
of the order parameter are developed simultaneously. 2) 
We construct a Ginzburg-Landau (GL) functional that 
contains a linear coupling between impurity spins and 
the d x 2_ y 2 + id xy order parameter. The time reversal 
violation is natural in this case as it allows the order 
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parameter to couple directly to the impurity spin. This 
coupling is possible only for d + id and not for d + is sym- 
metry of the order parameter. From the GL description 
we find that instability develops as a second order phase 
transition where both the out-of-plane magnetization and 
d xy component developed together and are proportional 
to each other. 3) This instability occurs at a certain con- 
centration range and is most likely suppressed at higher 
concentrations. 

Recent experiments reported the anomaly in the ther- 
mal conductivity in Bi2212 at low temperatures: the 
thermal conductivity of the Bi2212 with Ni impurities 
was observed to have a sharp reduction at T* = 200mK 
H . These experimental data indicate the possible super- 
conducting phase transition in the Bi2212 in the pres- 
ence of the magnetic impurities in a certain concentra- 
tion range. So far the transition has been seen only in 
the samples with magnetic impurities, e.g. Ni as op- 
posed to the nonmagnetic impurities such as Zn [0. It 
was reported that the feature in the thermal conduc- 
tivity is completely suppressed by applying the field of 
H ~ 300Gauss. The low field and the fact that feature 
disappears is consistent with the superconducting tran- 
sition into second phase. 

Our results might be relevant for the experimentally 
observed transition at T* in Bi2212 with Ni. 

We will construct a GL functional for the d-wave su- 
perconductor with magnetic impurities. Ni impurities is 
just one example of the magnetic impurities substituting 
Cu in Cu-0 planes. No other assumptions will be made 
about the nature of impurity but the existence of the non- 
compensated isotropic moment at low temperatures. We 
assume that superconductor is 2D with weak Josephson 
coupling between layers in Bi2212. The relevant fields to 
enter the GL functional are: 

• "to -the order parameter of the original d x 2_ y 2 com- 
ponent Ao = ^>o(cosk x a — cos k y a) , 

• $1- the order parameter of the d xy component Ai = 
*i sin k x a sin k y a 0, and 

• M- the magnetization of the impurity spins. 

Ao,i are the respective gap functions and we wrote 
them down explicitly in the momentum space, a is the 
lattice spacing. The impurity spin magnetization field 
M(r) = ^S(r) = £. ^S( r< )/(r - r<) with normal- 
ized to unity coarse-graining function /(r) which has fi- 
nite range , greater than the superconducting correlation 
length £ = Vf/Ao- This step is done to introduce a con- 
tinuum spin field in GL functional S. 

We consider a secondary transition \&o —> *o + ^i at 
T*. The relative phase of with respect to the phase of 
ty a is not fixed and will be determined by the free energy 
minimization. We will assume that the second transition, 
if at all, occurs at T* <C T c , where T c ~ 90K is the first 
transition temperature. Hence the order parameter 
which we can assume to be real, is robust and its free 
energy F(^q) can not be expanded in <S>q. 



Apart from coupling to the magnetic field B, which has 
to be determined selfconsistently, we assume the spin- 
orbit coupling between the impurity spins and orbital 
moment of the pair: 

H int = gL z M z (3) 

with some coupling constant g g|. Upon solving the 
microscopic selfconsistency equation one finds that this 
term generates the coupling between ^ and (see F int 
below) . 

Assuming expansion in powers of small M, \&i near 
second transition, the GL functional F = F sc + F magn + 
F int is: 

F sc = F(* 0) A) + + a 2 /4|*i| 4 + 

(3\{V i -2e/cA i )^ 1 \ 2 ,a 1 , 2 >0 

Fmagn = ^\M Z \ 2 + f\M z \ 4 + f |V;M Z | 2 + f£ 

F int = -£(*S*i - h.c)M z - (M z + - h.c.))B z (4) 

In writing the coupling between M z and ^g^i we nee d 
only the term which transforms as V^o and any function 
in the form ^o^Kl^ol 2 ) is allowed in the interaction term 
flO| . We will absorb $(|* | 2 ) into coupling b in F int , 
since is constant. The same absorption was done for 
the term with coupling v. should enter in Fi nt to sat- 
isfy the global U(l) symmetry: free energy should be in- 
variant under simultaneous \&o,i - * ^o,i exp(i#). F int is 
similar to linear coupling between magnetic field and the 
secondary order parameter, considered by Laughlin ||. 
Because of strong 2D anisotropy of Bi2212 we consider 
here 2D superconductor within the mean field approxi- 
mation. Weak interplanar coupling can be accounted for 
later. 

Below we consider Meissner phase, assuming that there 
is no net B and no mass currents in the bulk. This re- 
quires B = V T (M Z + ^(*o^i - h.c.)) = and can be 
satisfied for constant M z and L z . There is no coupling 
between L z and M z via magnetic field and there is only 
contact spin-orbit coupling. Only first term in Fi nt is 
nonzero therefore. 

All but F^t terms in the free energy Eq. (||) are positive 
and can not produce the instability of the original d x i_ y i 
state. Fi n t can be negative since it is linear in and M z 
and this term is crucial in producing second transition. 

In F mag n coefficients ai j2 are temperature dependent 
even at low temperatures. Direct exchange between Ni 
spins can be ignored in low concentration limit when the 
Ni-Ni distance is 20-30 A. The only relevant contribu- 
tion will be entropic term in F ma gn = E ma gnt - TS magn - 
In the presence of the ordered moments we get S magn = 
Sq — bnimpM 2 with some constant b. Hence a± — bTnimp 
in Fmagn- Rapid decrease of a\ with temperature might 
explain why second transition occurs at low tempera- 
tures. 
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Consider F sc . The second and third terms in F sc de- 
scribe be the cost of opening the fully gapped state with 
"Ji when interaction prefers to keep node, i.e. pure 
d x 2_ y 2 state. The change in free energy due to sec- 
ondary order parameter is given by the difference in en- 
ergy of quasiparticles before and after Hfi component is 
generated. We calculate the change in the energy of the 
superconductor subjected to the homogeneous external 
d xy source field: H xy = KEk^i( k ) c L c -k- ff ' where 
< K < 1 is the integration constant and Ck CT - elec- 
tron operator. Using standard result d K F sc — 1/ n{H xy ) 
we find increase of energy at Ai << T « A : 

SF SC = N Jo 1 dun J ^|Ai(0)| 2 In ( | Ao(0) |2 +K 2|i c l(0) |2 +T 2 )1 /2 
= lH*i| 2 + ^l*i| 4 

Here we take, based on the results of microscopic calcula- 
tion for the induced component, that is constant [ pT| . 
Nq is the Density of States at the Fermi energy and E c 
is the energy cutoff. 

The most interesting is the interaction term Fi nt . Few 
comments are in order here: 

1) The nature of the linear coupling between M z and 
Im^>}^i has simple physical meaning: for \E'o + *^'i state 
the orbital part of the pair wave function has a finite 
component of L z — 2 state: 

A ke _F5 = A + iAi = * cos 29 + i^x sin 29 

= ^ 1 exp{2iff) + (#o ~ *i) cos 29 (6) 

and interaction term describes the coupling between the 
orbital moment density L z oc Im^^i of the pair and 
the spin magnetization M z . One can show the existence 
of the linear coupling in microscopic theory by looking 
at the matrix element of Hi nt between and \&i states: 
Fi n t oc (*ocos26>| -ihd 9 M z \y i sin 29), L z = -ihd e , 9 is 
a planar angle. 

2) This form of interaction is possible only when \&i 
is the amplitude of the d xy component. 7m\f f g\f , i trans- 
forms as a pseudovector (magnetization) along z-axis. In- 
deed under time reversal T: ^0,1 ^$01 an d M x — > — M z 
so that Fi n t remains invariant. To check spatial inversion, 
consider it rotation around x-axis: x — > x, y — > —y, z —> 
—z. Under this operation M z changes sign, ^0 oc d x 2_ y 2 
does not change and ^1 oc d xy does. Hence the 7m$o$i 
changes sign and Fi n t remain invariant. One can check 
analogously other spatial transformations to see that in- 
deed Fi n t is a full scalar, d + is can not couple linearly 
to magnetic field M z . 

3) Linear coupling Im^^ \M Z enables the second 
phase transition, where both $1 and M z are developed 
below some critical temperature T*. 

4) The only orbital current allowed in the planar ge- 
ometry generate the z-axis magnetic moment which can 
couple only to M z . Since only Fi n t can be negative and 



drive the transition by coupling to M z , we can ignore 
other terms and consider free energy as a functional of 
M z only. 

Next, we will minimize the functional Eq.(Q) and find 
the solution for the ordered phase [ jl2| . Fix *5? a to be 
real positive and the relative phase of <3>i = |*i| exp(iv) 
without loss of generality. The set of equations is: 

■§j£ = aiM z - b sin ^0*1 = 
|£ = -&cos^ |* 1 |M 2 = 
^ = ai|*i| +a 2 |*i| 3 - bV M z smv = 

(7) 

From the second Eq.(0) it is obvious that v — ±ir/2 
and the phase will be determined by minimization of en- 

v = tt/2 sgn(bM z ) (8) 

This choice takes the maximum advantage of the L Z M Z 
coupling and minimization requires a complex order pa- 
rameter *o + i^i in the low temperature phase. T and P 
are violated spontaneously. Solving the remaining equa- 
tions we get: 

M z = £sini/tfo|¥i| 

I*iI 2 = ^(£*o-«i) = x(t;-t) 

5F= -^|*i| 4 ~ \T-T*\ 2 (9) 

This is the main result of this paper. Solution Eq.(|^) 
indicates that the transition is of the second order 
with the jump in the specific heat. We assumed that 
|\5 '$(T)\ 2 1 'ai(T) has a linear temperature slope near T*. 
It follows from the solution Eq(^) that: 

1) If coupling b is small enough b 2 < a\a\{T = 
0)/*q(T = 0) there is no nontrivial solution for For 
large enough b or small ai the solution near T* is de- 
termined by the temperature dependent ^q(T), ai(T): 
^q{T*) / a\{T*) — ai/b 2 . Strong temperature and con- 
centration dependence of a\ oc ni mp T might explain why 
the second transition occurs at low temperatures: the 
rapid decrease of a\ with T makes the condition for pos- 
itive solution in Eq(||) possible to satisfy at low T. 

2) As the function of impurity concentration two ef- 
fects occur simultaneously: condensate density l^ol 2 de- 
creases and a\ increases at increased impurity concen- 
tration. This effect will lead eventually to the disap- 
pearance of the transition when ot\ reaches critical value 
*q(T = 0) < a\a\/b 2 . Secondly, the suppression of *i 
due to increased impurity scattering will also lower T*. 
Quick suppression of the transition temperature T* with 
impurity concentration should be expected. 

3) The above calculation, which we have done in the ho- 
mogeneous approximation, is a variational proof of the 
instability of the original d x 2_ y 2 state toward the for- 
mation of the id xy component and spin magnetization 
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simultaneously. The order parameter and magnetization [1] 
in the real system do not have to be homogeneous. It is 
possible that there is inhomogeneous solution with lower 
energy. That is why the homogeneous solution should 
be considered as a variational approximation to the true ^ 
energy. 

4) Strong magnetic field parallel to the layers, H ^> 
H c i,ab ~ 1 Gauss in plane, will suppress the second 
phase. In the field Ni spins will be aligned in the lay- 
ers, linear coupling term on Hi nt will be zero and d xy 
component will vanish. This effect might explain the 
suppression of the second transition by magnetic field 
H c i ab *C H < H c i c ~ 300 Gauss , seen in experiment 

5) For real systems the interlayer coupling will lead to 
the true superconducting instability as opposed to the 2D 
mean field theory, considered here. Therefore details of 
interlayer interactions will be important in understand- 
ing true instability. It could be, for example, that real 
transition temperature is strongly determined by inter- 
layer coupling. [9] 

To test the proposed state following experiments can 
be done. The driving mechanism for the second phase 
clearly distinguishes between magnetic and nonmagnetic 
impurities, hence more experiments on Bi2212 with non- 
magnetic impurities will be helpful Because theory 
predicts the ordering of Ni moments below T* , one should 
be able to detect magnetization in uSR experiments. 
The increased superfluid density due to second compo- 
nent translates into the change in the penetration depth 
below T* which can be detected. These and other exper- 
iments will help to resolve if the proposed mechanism is 
correct. 

An interesting and related case is when Ni moments 
order intrinsically at some T m : a\ oc T — T m . The 
linear coupling, discussed here, with necessity will gener- 
ate d xy component at some temperature. It follows from 
the solution Eq. (||) with arbitrarily small a\ . First order 
transition is also possible. 

In conclusion, I presented the mechanism for a sec- 
ond order phase transition of original d-wave state into 
d + id state with spontaneously broken T and P. In the 
ordered phase both impurity spins M z and d xy compo- 
nent of the order parameter develop and are proportional [12] 
to each other. The low temperature phase develops mag- 
netic moment both due to magnetic impurities spin and 
because of the finite angular momentum of d + id state. 
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This work was supported by US DOE. 
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In the opposite case of the component open only in 
the narrow vicinity of the node, within width w oc ^i 
the angular integral in Eq.Q is restricted to small an- 
gles near the nodes. We find in this case that F sc oc 
l^i | 3 ln(_B c /max(^I'i, T)) which is similar, apart from log- 
arithmic factor, to the result obtained by Laughlin ||. 
We will focus on the ferromagnetically ordered case. The 
realistic configuration of impurity spins can be different 
and involve modulations, such as antiferromagnetic order 
within planes, alternating order within each plane, etc. 
These details depend on particular model one choses to 
address coupling between spins and conduction electrons. 
The goal here is to prove first that there is a variational 
state with total energy lower than pure d x 2_ y i state once 
spins are ordered. The lowest energy state might have in- 
homogeneous spin ordering and inhomogeneous supercon- 
ducting order parameter. 
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